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We discuss the phase diagram of the two-orbital model of half-doped manganites by calculating 
self-consistently the Jahn- Teller (JT) distortion patterns, charge, orbital and magnetic order at zero 
temperature. We analyse the instabilities of these phases caused by electron or hole doping away 
from half-doping, or by the application of a magnetic-field. For the CE insulating phase of half- 
doped manganites, in the intermediate JT coupling regime, we show that there is a competition 
between canting of spins (which promotes mobile carriers) and polaronic self-trapping of carriers 
by JT defects. This results in a marked particle-hole asymmetry, with canting winning only on 
the electron doped side of half-doping. We also show that the CE phase undergoes a first-order 
transition to a ferromagnetic metallic phase when a magnetic-field is applied, with abrupt changes 
in the lattice distortion patterns. We discuss the factors that govern the intriguingly small scale of 
the transition fields. We argue that the ferromagnetic metallic phases involved have two types of 
charge carriers, localised and band-like, leading to an effective two-fluid model. 



I. INTRODUCTION 

"Half-doped" manganites, corresponding to the gen- 
eral formula Rei_ x A x Mn03 with x = 1/2 where Re is 
a 3+ rare-earth metal ion and A a 2+ alkaline earth 
metal ion, eg., Lai/ 2 Cai/ 2 Mn03, have been the object of 
extensive experimental and theoretical studies for many 
years m& Here each Mn has an average valence of 3.5+ i.e., 
an average configuration of <i 3 ' 5 , corresponding to one 
Mn-e 9 electron for every two Mn sites hopping around 
amongst the two [(x 2 — y 2 ) and (3z 2 — r 2 )] e g orbitals 
on each Mn. The remaining three ti g electrons on each 
Mn are spin-aligned by strong correlations (Hund's rules) 
to form "core spins" with S — 3/2. Similarly to the end 
members, i.e., LaMnC>3 or CaMnOs, the half-doped com- 
pounds, thanks to their commensurate filling, are sim- 
pler in some ways than the doped manganites for gen- 
eral xi^ Nevertheless, they exhibit a very rich variety of 
properties as welli Specifically, as the system is cooled, 
there are phase transitions involving changes in magnetic, 
charge and orbital order, and, in some cases, metalicity. 
The details vary from material to material, depending 
systematically on the sizes of the "A site" ions of the 
perovskite structure. Generally, the lowest temperature 
phase is insulating, with simultaneous charge, orbital and 
CE type antiferromagnetic order (see below), and the 
charge/orbital order sets in first, at higher temperatures, 
compared to the antiferromagnetic order (e.g., for PrCa 
Tco/OO ~ 240K, whereas T N ~ 170122). The NdSr and 
PrSr systems show ferromagnetic metallic order at inter- 
mediate temperatures, but in the LaCa and PrCa sys- 
tems, the different phases obtained with increasing tem- 
perature continue to be insulating. Typically, the charge 
order and insulating behaviour at low temperatures per- 
sist on the "over-doped" side (x > 1/2), whereas the 
charge order disappears rather quickly on the "under- 



doped" side (x < 1/2), and is often accompanied by 
metalicity (except in the PrCa system, which stays in- 
sulating for all x). This asymmetry between "electron 
doping" and "hole doping" away from half-doping is a 
striking feature of the majority of the half-doped man- 
ganites. 

One simplifying feature of the half-doped manganites 
is that the low temperature phase is generally regarded 
as reasonably well characterised. Early neutron diffrac- 
tion work by Wollan and Koehler 4 suggested that the 
magnetic structure of La^Ca^MnOs can be viewed as 
a set of ferromagnetic zig-zag chains antiferromagneti- 
cally ordered relative to each other, with an 8-sublattice, 
(2-^/2 x 2-\/2) unit cell, and is referred to as the CE mag- 
netic order (Fig. HJ. The structure was qualitatively 
explained soon thereafter by Goodenough^ who pro- 
posed additionally that the phase also has a 2-sublattice, 
y/2 x \/2 charge order with alternating Mn 3+ and Mn 4+ 
sites, and a 4-sublattice, 2y/2 x y/2 striped orbital or- 
der, as indicated in Fig. Since then, CE order has 
been found in several other half-doped systems, such as 
Nd 1/2 Sr 1/2 Mn0 3 6 ' 7 or Nd 1/2 Ca 1/2 Mn0 5 8 , though some, 
such as Pri/ 2 Sri/ 2 Mn03)Sii show A-type antiferromag- 
netism, corresponding to [0, 0, 7r] spin order, i.e, ferro- 
magnetic planes of spins which are antiferromagnetically 
aligned in the z-direction). 

The presence of charge and orbital order is, however, 
harder to establish directly experimentally because of 
the lack of experimental probes that couple directly to 
them. Indeed, the perfect Mn 3+ /Mn 4+ charge order- 
ing proposed by Goodenough£ is currently regarded as 
controversia&ifiiiiiiS. X-ray diffraction data do suggest 
the presence of large Jahn- Teller (JT) distortions of the 
oxygen octahedra surrounding the Mn sitesii 3 - with two 
inequivalcnt Mn sites, of effective valence 3.5 + 6 and 
3.5 — 8, but 6 is not really known, and is unlikely to be 
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close to 0.5. In Pro.6Cao.4Mn03 (which is slightly under- 
doped though) , the charge or valence contrast seems fur- 
ther reduced, and it has been suggested on the basis 
of neutron diffraction data that the electron is rather 
shared by two Mn sites paired in dimer-like structures, 
referred to as "Zener Polarons", with (6 -C 0.5)^ It is 
not clear, however, whether this is due to the presence of 
additional electrons^ or a general feature of many half- 
doped manganites^ A recent work on Pro.sCao.sMnOa 
claims indeed to confirm the picture of the original CE 
state at precisely half-fillingji& An alternate, closely re- 
lated picture of the half-doped system is that of a bond- 
charge-density-wave, with no charge contrast of the Mn 
ions, substantial hole occupancy on the oxygen ions on 
the chains, and alternating "0 2 ~/0~" order. 17 There 
have been X-ray studies on orbital order and correlations 
as well as charge and magnetic order using a variety of 
methods such as soft x-ray resonant diffraction^ coher- 
ent x-ray scattering^ which explore the spatial extent 
of orbital order, in particular. The resonant scattering 
experiments in Pro.6Cao.4Mn03 conclude that the charge 
disproportionation is less than complete, that there is or- 
bital mixing, and that therefore the simple Goodenough 
model is not right. On the other hand, when holes are in 
excess (x > 0.5), it has been suggested that charge order 
persists but becomes incommensuratei22*2I 

A closely related family of manganites carefully stud- 
ied recently is Ao.sA'o.sMnOa where A is a rare earth (Y, 
Tb, Sm, Nd, Pr, La in order of ion size) and A' is Ba^ 
The phases have been studied as a function of A- A' site 
order/disorder. When there is order, the low tempera- 
ture phase is charge-ordered (CO) and orbitally ordered 
(00) for ion size from Y to Nd. Beyond Nd, up to La, 
the phase is a ferromagnetic metal. However, when A-A' 
sites are disordered, there is no CO phase at all, but only 
a spin glass (SG) phase, from Y to Sm, after which the 
ground state is ferromagnetic (FM). This means that the 
CO/FM and CO/SG competition depends on ion size as 
well on A site ordering in the perovskite ABO3 structure. 

Another interesting and intriguing property of the 
CE charge-ordered (CE-CO) phase is the magnetic- 
field-induced insulator-metal transition first discovered 
in (Nd,Sm) 1/2 Sr 1/2 Mn0 3f 22i2i and later shown to be 

ubiquitous^. Though insulating at zero field, these ma- 
terials become metallic by the application of magnetic- 
fields of the order of 10 - 40 Tesla via sharp, first- 
order, resistive transitions^ The magnetic-field ener- 
gies involved are much smaller than the thermal ener- 
gies (of order 200 K) needed to destroy the charge or- 
der, and orders of magnitude smaller than the charge 
gap of 0.2 - 0.3 eV, as observed as a function of field 
by STM in Ndo. 5 Sro.5Mn0 3 i2£ This can be viewed as a 
different manifestation of the colossal magneto-resistance 
seen at the metal-insulator transition of doped mangan- 
ites for x ~ 0.25ji and the microscopic understanding 
of the above features poses similar difficult theoretical 
challenges i£ 

Theory. A variety of models and mechanisms have 




FIG. 1: (color online). A depiction of the CE charge-ordered 
antiferromagnetic phase (CE-CO). The "bridge sites" (1,3), 
at the centers of the arms of the zig-zag chains which are 
ferromagnetically ordered, have alternate occupancies of the 
3x 2 — r 2 and 3y 2 — r 2 orbitals. At the "corner sites" (2,4), 
there is no orbital order, unless JT interactions are present. 
5 is the charge disproportionation between the occupancies of 
the corner sites and the bridge sites. 



been examined in the context of half-doped manganites 
as wellii The simplest model has mobile electrons mov- 
ing amongst non- degenerate orbitals, coupled to the Mn 
tig core spins by a large Hund's rule (double) exchange 
coupling Jh- The latter promotes ferromagnetism, but 
competes with a direct antiferromagnetic coupling Jaf 
between the core spins. Even in this simple model, the 
ferromagnetic or CE types of order are stabilised depend- 
ing upon the strength of Jaf Van den Brink et at 
considered a more realistic model with the two types 
of e g orbitals of Mn, and found that the CE phase is 
orbitally-ordered: the " bridge sites" of the zig-zag chains 
have alternating preferred occupancy of (3x 2 — r 2 ) and 
(3y 2 — r 2 ) orbitals^ as indicated in Fig. ^ They also 
showed that a charge contrast 8 (not bigger than 0.2) 
can be generated by including on-site Coulomb inter- 
action. This is because the "corner sites" turn out to 
have equal occupancy of (x 2 — y 2 ) and (3z 2 — r 2 ) or- 
bitals, and that costs Coulomb energy. To reduce this, 
the system adopts a preferred occupancy of the bridge 
sites which are orbitally ordered. In later work, nearest 
neighbour Coulomb interactions were also includedi^iSl 
However, the charge order due to long-range Coulomb 
interactions is generally of the Wigner-type with wave 
vector Q — (tt, tt, tt), contrary to the (tt, tt, 0) order, with 
charge stacking along the z-direction, suggested by exper- 
iment. To stabilise the (tt, tt, 0) order in a wider regime 
of parameters, JT interactions between the Mn ions and 
their surrounding oxygen octahedra, which are supposed 
to be quite large^ have to be included^ The conse- 
quent JT distortions further lower the energy of the CE 
phase because of the already present (2>x 2 — r 2 ) / (3y 2 — r 2 ) 
orbital order. Classical Monte Carlo simulations includ- 
ing static JT distortions on small clusters as well as self- 
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consistent mean field treatments of models including JT 
and Coulomb interactions 3 - 1 suggest that the CE charge 
stacked state has the lowest energy in an intermediate 
range of Jaf, unless the nearest neighbour coulomb in- 
teraction V becomes much too large. 

However, to our knowledge very few of these studies 
have addressed the other issues, such as the magnetic- 
field-induced insulator metal transition, and the clcctron- 
hole asym metry. The first issue was tackled in refs. 
|2fil27l32| by assuming model parameters very close to 
the phase boundary between the ferromagnetic and CE 
states. The resulting small energy difference between the 
two phases can then be overcome by an arbitrary small 
magnetic-field. But it is hard to justify why the sys- 
tem parameters should be so finely tuned for so many 
systems. As regards the second issue, band structure 
arguments r& and treatments including JT distortions on 
small clusters^ necessarily lead to metallic phases upon 
addition of electrons or holes, in contrast to experiments. 

Recently, a theory for doped manganites has been 
proposed^ where it is argued that due to strong JT in- 
teractions the e g electrons dynamically reorganise them- 
selves into two types. The majority of the electrons (la- 
belled £) become localised polarons, trapped by large lo- 
cal JT distortions; and a minority of them (labelled b) can 
remain mobile and non-polaronic. Still virtual adiabatic 
transitions to empty neighbouring sites induces a ferro- 
magnetic exchange referred as virtual double exchange. 33 
The resulting Falicov-Kimball like, I — b model Hamilto- 
nian treated in a simple dynamical mean-field treatment 
in the framework of an "orbital liquid" description, gave 
a good account of the properties of doped manganites— 

In this paper, we propose an extension of the above 
theory to the half doped case, which has to include the 
possibilities for orbital, charge, and antiferromagnetic or- 
der. We obtain pointers to this by studying the proper- 
ties of electronic excitations coupled with JT defects in 
the lattice distortion pattern. We find that such a study 
suggests the incipient instabilities of the CE phase in- 
dicative of the doping and magnetic-field induced phase 
transitions seen experimentally, as well as the presence 
of localised and mobile carriers. The localised states on 
the defects which we obtain are different from the micro- 
ferrons suggested at small x around a dopant^ as they 
are self-generated and could exist even in the absence of 
chemical disorder. In principle, the JT defects we are 
considering could be mobile on a longer time-scale, al- 
though disorder may indeed pin them down. 

More specifically, in this paper, we first determine the 
zero temperature phase diagram of the 3d two-orbital 
model of half-doped manganites for periodic phases in 
the thermodynamic limit, including JT distortions, but 
ignoring Coulomb interactions. We do this by minimising 
the energy assuming a periodic unit-cell of eight sites^i 
inside which static JT distortions and core spin directions 
are allowed to be arbitrary. This allows us to determine 
them self-consistently without using finite-size clusters, 
thereby extending and reinforcing earlier work— iSl In 



particular we obtain analytic results for the phase bound- 
aries at strong JT coupling. 

Next, we show that the periodic ferromagnetic phase 
obtained at small Jaf by the method discussed above can 
become unstable with respect to a phase with two types 
(I — b) of electrons when the JT coupling is lowered. We 
show indeed that it becomes energetically favourable to 
create single site JT defects, i.e. release the distortion 
on a finite number of sites that were previously distorted 
and promote previously trapped electrons onto a mobile 
band, thus suggesting a metallic phase. The exact nature 
of the phase can not be figured out by such an instabil- 
ity analysis. Nonetheless, it suggests an effective (I — b) 
Hamiltonian with orbital degrees of freedom explicitly 
included. 

The observed phases at half-doping, such as the CE 
phase, are antiferromagnetic, corresponding to appropri- 
ately larger values of Jaf- But they show transitions to 
ferromagnetic metallic phases in an external magnetic- 
field or when electrons are added. To understand such 
transitions, in addition to considering changes in the JT 
distortions, canting of spins is important. 

Canted phases are expected to appear not only in a 
magnetic-field, but also upon doping with carriers (and 
irrespective of their nature), following the original ar- 
gument by de Gennes— Here we show, however, that 
canted metallic phases appear only when electrons (and 
not holes) are added, because of the underlying asymme- 
try of the density of states at half-doping. When we allow 
for JT distortions, we find a competition with a disor- 
dered phase where the added electrons are trapped by JT 
distortions, the latter phase winning only at small elec- 
tron concentration. On the hole-doped side, added holes 
are simply trapped by the lattice distortions and the sys- 
tem remains insulating. Thus our work provides an ex- 
planation for the particle-hole asymmetry near x ~ 1/2, 
at intermediate JT couplings which we argue are relevant 
for the majority of the manganites. 

Similarly, we find that there is a strong interplay be- 
tween turning on a magnetic-field at half doping and the 
JT distortion pattern. This is consistent with x-ray mea- 
surements in Lai/2Cai/2MnC>3 in a field— i£i Starting 
from the distorted CE phase, we find in addition an in- 
stability of the high-field ferromagnetic phase to the for- 
mation of JT defects. This suggests that the high-field 
phase seen in experiments may need a two-fluid descrip- 
tion. 

Interestingly, a very similar two-carrier hypothesis was 
proposed based on phenomenological grounds to under- 
stand the resistivity of Lai^Ca^MnOs (x ~ 1/2). 36 
More recently, a particle-hole asymmetric Ginzburg- 
Landau theory was proposed to explain^ the incommen- 
surate charge order— seen for x > 0.5. We believe that 
our theory provides the microscopic basis for these facts 
both. 

The rest of this paper is organised as follows. In sec- 
tion [HI we discuss the phase diagram of the half-doped 
manganites restricted to periodic ground states with the 
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most general 8-sublattice structure. We give in particular 
an analytic strong-coupling description 1(11 D(l . In section 
IIIII we study the instabilities of some of these phases : in- 
stability of the strong JT coupling ferromagnetic phase, 
which defines a new phase (|III Af) ; instability upon dop- 
ing to the canted phases (fill B 1(1 or to self-trapping of 
added carriers (fill B 2(1 , and the competition between the 
two l(IIIB3() . The effect of the magnetic- field on the CE 
phase is studied in section ITlI CI where we also discuss 
the nature of the high-field ferromagnetic phase. In sec- 
tion IIVI we summarise and discuss the possibilities for 
testing these ideas experimentally. A short account of 
some of these results has been presented elsewhere^ 

II. OPTIMISED PERIODIC PHASES AND 
PHASE DIAGRAM FOR HALF DOPED 
MANGANITES 

A. Model Hamiltonian 

Our discussions are based on the following Hamiltonian 
for the manganites 

H[{S la ,Q ia ,Q la }} = - ^2 Kbi](SiaiS b)c\ aa C jb {3 
ijct{3ab 

+ JAFSta.Sjfe — gjlB ^] H ■ Sj a 

<ijab> ia 

+ \ K ^Qia -dYl QiaT a 0(Qia)4 aa Ci a (1) 
ia iaa/3 

where c\ aa creates an e g electron (in a low-energy- 
projected Wannier orbital with e g symmetry 41 ^, on the 
sublattice site a (Mn site) of the unit-cell i in the 3d 
cubic lattice, and in the orbital state a = 1,2, with 
1 = d x 2_ y 2,2 = d^ z 2_ r 2. The original lattice is de- 
composed into eight sublattices, labelled with There 
are N sites and cN — (1 — x)N electrons (when x = 1/2 
the number of electrons is denoted No = N/2). The first 
term is the kinetic energy of the electrons. The hop- 
ping parameters are taken to be of the usual Anderson- 
Hasegawa form 42 which takes care of the Hund's cou- 
pling, Jh J2i Sia ■ Si a in the limit of large Jn/t, with 
Sia, the S — 3/2 core spin formed from the Mn i 2g elec- 
trons being approximated as a classical spin. As a conse- 
quence only the electrons with spin projections parallel 
to the core spins are present, and their hopping ampli- 
tudes are functions of the polar angles of the core spins 
given by 424 s . 

lap ,„ s n _ ,a/3 
L abij \&ia,&jb) — ^ abi j * 

cos — cos — — h sin — sin— e w '° ^'(2) 
2 2 2 2 ) ' 

Here f"j^- is the usual, anisotropic and symmetry de- 
termined, hopping amplitude 31 between the e g orbitals 



a and j3 at the two nearest neighbour sites (i, a) and 
(J, b) respectively, arising from their hybridisation with 
the O — p a orbitals (with 4i/3 being the hopping be- 
tween (3z 2 — r 2 ) orbitals in the z-direction)2i. The sec- 
ond term Jaf is the antiferromagnetic coupling of the 
t2 g core spins that comes from standard superexchangc 
processes^ It can be roughly estimated from the Neel 
temperature of a system with only t^g core spins, such 
as CaMn03, although the structure of the half-doped 
system is not exactly the same. The third term is the 
Zeeman energy where H is the external magnetic-field. 
The last two terms include the Jahn- Teller (JT) phonons 
and their coupling to the e g electrons. We neglect the 
Pf a /2M ia terms since huj <C t (where ojq is the typical 
phonon frequency), but include their effects heuristically 
when we argue that JT defects lead to polaron formation. 
Qia and Qi a are, respectively, the amplitude (measured 
in units of the typical JT distortions in these materials) 
and the angle of the JT distortion at the site (i, a). The 
coupling matrix is given by: 

/cosG sine \ 

v 1 \ sm e — cos e J w 

K is the lattice stiffness of a simplified non-cooperative 
model where distortions on neighbouring sites are not 
coupled. More detailed and realistic models would in- 
clude cooperative JT couplings and couplin g to breathing 
modes such as in the lattice model of Ref. |46j . 

We have neglected the on-site Coulomb interaction, 
U J2i n iaaniaa between different orbital states. Although 
it is an important interaction in the problem, we can not 
treat it using the methods used in this paper except in 
a mean-field approximation. However, it is effectively 
taken into account when orbital order occurs. We com- 
ment on the effects of its inclusion at appropriate places 
in the paper. When a local JT distortion occurs on a 
site, the degeneracy of the e g orbitals is lifted. If only 
one electron is present, there is a gain by occupying the 
lowest energy level. If a second electron is added, how- 
ever, it has to occupy the higher energy level because 
of the strong Hund's coupling. There is a compensation 
and the energy gain vanishes. JT distortions therefore 
suppress double occupancy of sites, mimicking the effect 
of U. When g <C t, the distortions are small or zero and 
it is important to explicitly include U, which does play 
a role. For instance, it induces a charge-ordering in the 
CE phased just as a finite g/t does^i When g 3> t (see 
section Til Dll . JT distorted phases appear naturally, and 
the inclusion of U is not crucial. Similarly, the inclusion 
of the term U n^nii IS unimportant (completely irrel- 
evant when Jh — > 00) because the large Hund's coupling 
prevents double occupancy of this type. We have also 
neglected the longer range coulomb interactions as they 
are expected to be weak because of the large dielectric 
constant of the manganites, and we do not consider is- 
sues (such as macroscopic phase separation) which are 
sensitive to their presence. 
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Note also that regarding the direct coupling of the ti g 
spins, we restrict ourselves to a pure Heisenberg superex- 
change coupling. This may not be absolutely accurate 
for S = 3/2 spins. Further couplings, such as single-ion 
anisotropics, are certainly present in the real materials, 
but are not of crucial importance for the issues we focus 
on in this paper. We therefore restrict ourselves to the 
Hamiltonian of eq. Q. 



B. Method 

The Hamiltonian of eq. (JIJ represents mobile electrons 
coupled to local classical degrees of freedom that act like 
annealed disorder. The probability weight of a partic- 
ular configuration of the classical variables is given by 
exp[— F e i/(kBT)] where F e i is the electronic free energy 
in the presence of that configuration; their distribution 
thus has to be determined self-consistently. At zero tem- 
perature, it is reasonable to assume that the classical de- 
grees of freedom are frozen and have well-defined values. 
On a finite lattice these can, in principle, be determined 
as follows. One can diagonalise the Hamiltonian exactly 
for a given (arbitrary) configuration of lattice distortions 
(Qi,®i) and polar angles of the spins 9i (for simplicity 
we are ignoring the azimuthal angles of the spins fa). For 
each configuration, one can thus find the electronic en- 
ergy levels, fill up the states up to the Fermi energy, and 
calculate the total energy. To obtain the ground state of 
the system, one then needs to minimise this energy with 
respect to all the possible configurations of classical vari- 
ables. Such a procedure can be implemented, for exam- 
ple, using a Monte-Carlo technique-^ for a finite lattice, 
but becomes a more and more difficult task as the num- 
ber of lattice sites, and hence the number of variables, 
increases. 

Since our aim is to explore the physics of the exper- 
imentally observed CE state, we adopt a simpler ap- 
proach. We assume a 8-sublattice periodic structure that 
is compatible with the periodicity of the CE state, which 
permits us to tackle the problem in a lower dimensional 
space of the classical variables. [Needless to say, this 
rules out the possibility of incommensurate (with respect 
to the assumed eight sublattice structure) or inhomoge- 
neous phases.] We have implemented a simulated an- 
nealing routine to minimise the total energy with respect 
to (essentially all possible) distortions and spin angles 
on the 8-sublattices. Thus our approach differs from 
and is complementary to earlier numerical approaches 
which have considered small clusters and done a full 
classical Monte-Carlo simulation for the spin and lattice 
variables^ For us, the only limitation is the number of 
sublattices, which we fix to be eight; the system size is 
vastly larger (we are treating the real 3d case), and prac- 
tically in the thermodynamic limit. [The system size, 
i.e, the total number of sites, is 8 times the number of 
periodic repetitions of blocks of 8 sites (the spin and dis- 
tortions being the same in all the blocks), and we have 



done calculations using up to 6912 blocks.] The compu- 
tational effort of such an approach compared with that 
of Ref. |u] is, on one hand, much reduced because we do 
not have to equilibrate a large number of variables. On 
the other hand, the calculation of the energy of each con- 
figuration takes more time because we sum over a large 
number of k-points (equal to the number of blocks) in 
the Brillouin zone corresponding to the periodic struc- 
ture. We have carefully studied the finite-size effects on 
these Brillouin zone sums. The error on the total en- 
ergy due to the truncation of the sums is of the order of 
10~ 2 t where t is the typical energy scale of the problem. 
The thermodynamic limit is therefore reached within this 
accuracy, i.e., for all the phases compatible with the sub- 
lattice structure we expect that our results are within 1% 
of the thermodynamic limit results. 



C. Results 
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FIG. 2: Phase diagram of the 3D two-orbital model (T = 0, 
x = 0.5, K/t = 10). FM (resp. FM d ): ferromagnetic metallic 
phase with no distortions (resp. small uniform distortions). 
FI-CO (resp. FM-CO): charge- ordered ferromagnetic insulat- 
ing (resp. metallic) phase with distortions that favour occu- 
pancy of the x 2 — y 2 orbitals (Fig. |3J. A^: ferromagnetic 
planes AF aligned with uniform distortions. A-CO: A with 
charge order. CE-CO: Ferromagnetic zig-zag chains AF or- 
dered, charge and orbital ordered (3a; 2 — r 2 /3y 2 — r 2 ) [Fig. 
0. G-CO: Neel AF phase with charge-order. Inc.: Possible 
incommensurate states that interpolate between CE and G. 
Dotted dashed lines come from analytical expressions derived 
in the strong-coupling limit (section III Dll . Solid (dashed) 
lines show first-order (second-order) phase transitions. 

Figure [3 shows the phase diagram as a function of 
the JT coupling, g/t, and the antiferromagnetic cou- 
pling, JafS 2 /t, at zero temperature. [We choose units 
such that the JT distortions are dimensionless, whence 
K and g both have dimensions of energy, which we spec- 
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FIG. 3: (color online). A depiction of the ferromagnetic in- 
sulating (resp. metallic) charge-ordered phase (FI-CO [resp. 
FM-CO]) stable at strong JT coupling and small antiferro- 
magnetic coupling (see fig. |5J. Alternate sites have charge 
disproportionation S (given in Fig. |S]as a function of g/t). 
The sites with higher occupancies also have strong JT dis- 
tortions (see fig. 0J, of such orientation as to promote the 
occupancy only of (x 2 — y 2 ) on these sites, leading to orbital 
order as well. Note that the lattice is rotated by 45 degrees 
with respect to fig. Q 



ify in units of t. We use a fixed K/t = 10. To com- 
pare with earlier work, the JT energy is then Ejx/t = 
{g/t) 2 /(2K/t) = {g/t) 2 /20.} We basically find the same 
phases that were found before either by comparing the 
energies of selected phases at g = 0(22*22^ or by Monte- 
Carlo simulations at finite <?j2i except that now we have 
provided confirmation that they are indeed the optimal 
8-sublattice structures in the thermodynamic limit. 

We now describe the different phases shown in Fig. |3 
including the amplitudes of the JT distortions in them 
and the corresponding electronic properties. 

For small values of Jaf an undistorted metallic phase 
with 3-d ferromagnetic order (FM) is stable up to a crit- 
ical value of g/t ~ 5. Above this threshold, there is 
a ferromagnetic phase with very small uniform distor- 
tions (Fig. HJ, noted FM^. There is also a narrow region 
(5.6 < g/t < 5.9) where the solution displays many in- 
equivalent sites. As discussed later, we believe that this 
is consistent with the instability that we find in section 
ITTT1 For g > 5.9, the stable phase is the FM-CO fol- 
lowed by the FI-CO. These phases have the structure 
depicted in Fig. [21 i.e., a layered structure with large 
JT distortions on two sites out of four in a checker-board 
pattern in each layer, and oriented in such a way as to 
favour the (x 2 — y 2 ) orbital on the strongly distorted sites 
(Fig. 0). There is a charge disproportionation 8 that is 
given in Fig. [3J The structure is metallic (FM-CO) up 
to g/t ~ 6.3, and is insulating (FI-CO) for larger g/t, 
as is clear from Fig. |SJ This structure has been found 
before^ and is known to compete with a similar struc- 
ture which prefers a (3a; 2 — r 2 )/(3y 2 — r 2 ) orbital order, 
when strong anharmonic and cooperative JT couplings 
are taken into account^ 

The A phases, which are more stable at larger Jaf (the 
larger the g/t, the smaller the coupling Jaf required for 
the transition) are similar to the ferromagnetic phases we 



have just described except that successive layers are now 
antifcrromagnetically ordered. For g/t up to ~ 5.1, the 
phase noted A^ is uniformly distorted with a distortion 
amplitude given in Fig. It is metallic in this regime 
(see the charge gap in Fig. |HJ). For larger values of g/t, 
the A phase become charge-ordered and insulating, as in 
case of the FI-CO phase (Fig. E). 

The CE-CO phase, the CE phase with charge stacked 
order and orbital order (Fig. QJ, is the stablest over a 
wide range of parameters for intermediate Jaf, as is clear 
from (Fig. As pointed out in Refs. |28u3ll |. there is 
orbital ordering even at g — 0, but no charge ordering; 
the "bridge sites" having an average occupancy of 0.5, 
but only of (3a; 2 — r 2 ) and (2>y 2 — r 2 ) orbitals alternately 
(Fig. nj. The "corner sites" on the other hand, have 
equal occupancy (0.25 each) of {x 2 — y 2 ) and (3z 2 — r 2 ) 
orbitals. The sites are undistorted at g = 0, but get dis- 
torted as soon as g > (Fig. bottom- left panel). The 
distortions on the bridge sites are the largest, and are ori- 
ented in such a way as to further stabilise the alternating 
occupancy of the (3a; 2 — r 2 ) and (3y 2 — r 2 ) orbitals that 
already exists at g = 0, since distortions that precisely 
favour this alternation lower the energy of the system. 
In addition, small distortions get generated also on the 
corner sites that favour the (a; 2 — y 2 ) orbital (Fig. 0J. 
As a further consequence, a charge disproportionation 6 
between the bridge and corner sites develops, favouring 
a higher occupancy of the former. The variation of 5 
with g/t is shown in Fig. [SJ The system is an insulator 
whatever the charge disproportionation, as shown by the 
finite charge gap in Fig. 

For strong JT coupling (g/t 3> 1), the CE-CO phase 
is degenerate energetically with the C-CO phase which 
consists of straight ferromagnetic chains antiferromagnet- 
ically ordered with respect to each other. The charge or- 
der is accompanied by orbital order of the (3z 2 — r 2 ) type 
if the chains are oriented along the z-direction. This de- 
generacy is discussed in explicit detail in section Hi Dl In 
this limit, it is easy to show (see section Hi Dfl that the G 
phase (completely 3d-AF phase with localised electrons) 
is always the stablest for large values of Jaf- For smaller 
values of g/t (and large Jaf) we find solutions with many 
inequivalent (canted) sites, suggesting that the transition 
from the CE phase to the G phase in this regime might 
proceed via intermediate states that are incommensurate 
relative to the periodicity of the unit-cell we have consid- 
ered (Fig. 

We consider next the strong JT coupling limit, g/t 3> 
1, whence we can calculate the energetics of the phases 
and the phase boundaries discussed above analytically. 



D. Localised description for g/t > 1 

At large g/t, it is energetically favourable for all the 
electrons in the system to be self-trapped by local lattice 
distortions since the JT energy gain is large. We hence 
start with Wannier-type wave-functions with electrons 
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FIG. 4: Amplitudes of distortions of the four inequivalent sites 
of the unit-cell as function of the coupling parameter g/t. The 
different panels represent the various phases found previously: 
F,A,CE,G. The last three phases do not exist for all values 
of g/t; the curves are then obtained by fixing the magnetic 
structure and optimising with respect to the distortions. For 
the G-CO phase, the distortions are exactly given by g/K 
since the electrons are completely localised. 




g/t 

FIG. 5: Charge disproportionation defined by the valence of 
the two inequivalent Mn ions, Mn 3,5+i and Mn 3 ' 5 ~* (see Figs. 
0and[3J, as a function of g/t (or E.jr/t [inset] defined by 
IEjt = gQmax, where Q ma x is the distortion of the site with 
the largest distortion) for the CE, F and A type-phases. 



fully localised on strongly distorted sites. The local en- 
ergy per electron is the sum of the elastic energy \KQ 2 
and the electronic energy gain —gQ which is minimal at 
Q = g/K with a net energy gain of Ejt — g 2 /(2K). 
In the limit of large g/t, there is a large degeneracy be- 
cause electrons can be trapped on any site and in any 




FIG. 6: Charge gap vs. g/t for the CE, F and A type- phases. 
The CE phase is always insulating, while the A and F phases 
are insulating beyond g c A/t = 5.1, and g C F/t ~ 6.3. The 
dotted line is the gap obtained with non-optimised distortions 
(eq. JITJ in Elig . 



orbital state, as long as the orbital state correlates with 
the orientation 8, of the JT distortion as 



l*(e«)) 



cos— 14 



(4) 



(for instance, 0^ = — tc/3 for 3x 2 — r 2 , and 8, = ir/3 for 
3y 2 -r 2 ). 

This degeneracy is lifted at second-order in pertur- 
bation theory in the kinetic energy of the electrons. 
Consider an electron localised on a site with an empty 
neighbouring site and with the corresponding core spins 
aligned. Then, in the adiabatic limit (t ^> Tiujq where lvq 
is the frequency of the JT phonons) appropriate here, it 
can hop virtually onto any of the two orbital states of that 
site and back, without giving the lattice distortions time 
to relax (the relevant energy denominator being 2Ejt), 
and hence lower its energy. This energy lowering is less 
if the core spins are misaligned whence the hopping am- 
plitude is reduced (and even fully suppressed in case of 
anti-alignment). It is also less if the neighbouring site 
is occupied, whence the energy denominator is larger, 
equal to 4Ejt (4Sjt + U in the presence of U, so that 
the process gets suppressed altogether for large U). Such 
a process hence gives rise to a effective double exchange 
term in the Hamiltonian as pointed out in ref. |3^| and 
labelled virtual double exchange. The dominant term in 
the effective Hamiltonian is then4&: 



ft = -E., T J2 n ^+J2 J AFSi-Sj ~ SMb H S < 



J 



jE^' S ^ + M 1 - »J)C?J + (i " j)] (5) 
<i,j> 

Here Ejt = g 2 /(2K), m is the electronic occupancy 
on site i, J — t 2 /{2E JT S 2 ), i = 4t/3. C itj = 



8 



cos[(9i + %)/2] with f 



i,i-\-x 



^ i.i-\-y ~ ^i-\-y.i 



-7I-/3 and ^ iyi+z = 



-tt/3, 



The effective Hamiltonian of eq. is a classical spin- 
charge-orbital model with no quantum fluctuations. If 
the charges are assigned specific positions (so as to min- 
imise the energy), the model reduces to a spin-orbital 
model. It is different from the spin-orbital model pro- 
posed for undoped LaMnC>3 obtained by projecting out 
double occupancies^ because double occupancy is irrel- 
evant in the limit being explored here. The orbital (and 
JT distortion orientation) variables on neighbouring sites 
are not directly coupled in this model (as C^- involves 
only one orbital angle 6.;). Such a coupling would arise 
if one takes into account a direct coupling between JT 
distortions on neighbouring sites, as in the cooperative 
JT model. Nevertheless, even in this simplified approach, 
the virtual double exchange lifts the degeneracy between 
different orientations of the JT distortions and the cor- 
responding orbital degeneracy. In addition, it clearly 
favours ferromagnetic bonds and charge disproportion- 
ation. 

We can now understand the strong-coupling limit of 
the phase diagram (Fig. EJ) , and furthermore even in the 
presence of a magnetic field, by estimating the energies 
of the various phases using equation (0. For the fully 
charge disproportionated phases the energies per site are 
obtained by minimising TL with respect to the canting an- 
gle (at finite fields) and orbital angle (the latter depends 
on the field for the CE phase, and only the leading term 
in H — > is given): 
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At zero field, there is only one free parameter, Jaf I J, 
which determines the relative energies. At Jaf = 0, the 
FI phase (which is orbitally disordered in this limit) is de- 
generate with the A- type phase (with (x 2 —y 2 ) orbital or- 
der), but the latter is favoured as soon as Jaf > 0. There 
is a succession of first-order phase transitions as Jaf/ J 
is increased, first to the CE-CO phase (degenerate with 
the C-CO phase) at Jaf I J =1/4 and then to the G-CO 
phase (see Fig. 0. In terms of the original variables, the 
first transition at 1/4 is located at JafS 2 /t = AtK/ (9g 2 ); 
the second transition is at JAFS 2 /t = 8tK/(9g 2 ). The 
phase boundaries given by these equations in the strong- 
coupling regime are displayed as dash-dotted lines in the 
phase diagram (Fig. |2J) and are in good agreement with 
the phase boundaries obtained numerically. 



Similarly we can compare the energies of the different 
phases as a function of the magnetic-field and draw the 
corresponding phase diagram in the (Jaf/ J, Eh^bH/JS) 
plane (Fig. 0). The degeneracy between the C and the 
CE phases is lifted and the C phase wins at finite fields. 
This is because the fourth term of (jHJ favours Wigner- 
crystal type of ordering. For the C phase, for instance, 
the critical field is g/isHc = 8JafS — JS. 
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FIG. 7: Phase diagram when g/t > 1 at T = (x = 0.5). 
The phases are the same as in Fig. JS 2 = P/2Ejt. The 
CE-CO and C-CO phases are degenerate at zero-field, but the 
latter wins at finite fields. 

As noted earlier, in the limit of large g/t much of 
the physics is insensitive to the inclusion in the model 
of the on-site Coulomb interaction between different or- 
bital states, Un ia n i& . The total energy is hardly affected 
since double occupancy is much reduced. We emphasise 
again that this is contrary to what happens in the other 
limit j/i< 1 where the electron density is uniform. In 
the latter case, it is known that U by itself will induce 
charge-ordering in the CE phased at least if the latter 
is not destabilised by other phases^i In the CE phase 
at large g/t, U will slightly modify the charge contrast 
by pushing the electrons further off the corner sites. We 
have performed a self-consistent Hartree-Fock calculation 
to confirm this point. At g = 0, the calculation gives the 
same results as the slave-boson approach^ At g/t = 7 
and for the optimised lattice distortions, the charge con- 
trast increases, with respect to U — 0, by a very small 
amount of order 0.05 for U/t even as large as 25. There- 
fore, we conclude that the effect of U is small and does 
not change the nature of the insulating phases in the limit 
of large g/t. 

We next consider the interesting question as to what 
the appropriate low energy effective Hamiltonian replac- 
ing © is when t/EjT becomes sufficiently large that per- 
turbation theory in t/Ejx and the Hamiltonian JSJ are 
not valid anymore. We argue below, by studying the 
excitations and instabilities of the original model JQl in 
the ferromagnetic phase, that the effective model that 
replaces (0 when t/Ejx gets larger takes a similar form 
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except that mobile electrons have to be included. 

III. INSTABILITIES OF THE PERIODIC 
PHASES 

We have discussed in the previous section the vari- 
ous phases stable in the thermodynamic limit that are 
periodic with a 8-sublattice unit-cell. We will discuss 
in this section several instabilities that point to non- 
periodic phases, at half-doping (fill Afl and also upon dop- 
ing IjHI B|) or addition of a magnetic-field IjHI C|) . We will 
show that the ferromagnetic insulating phase (FI-CO) is 
in fact unstable when g /t is lowered below a critical value 
g c /t (~ 6.8. This instability occurs before any of the tran- 
sitions discussed above (at g c /t ~ 6.3 and 5.9) take place. 

For this purpose, we study the energetics of defects 
in the lattice distortion pattern of the periodic phases. 
We consider particle and hole excitations accompanied 
by single site JT defects. We consider both types of de- 
fects, one where we add a distortion on a site that was 
previously undistorted, and the other where we remove 
the distortion of a site that was distorted. Without the 
lattice distortion defect, the lowest energies of the par- 
ticle or hole excitations accessible are the appropriate 
gaps determined by the band-structures discussed in sec- 
tion HIO (Fig. El. The defect modifies locally the JT 
energy level, and hence constitutes a scattering potential 
for the particle and hole excitations. The problem lacks 
lattice translation invariance, and we have solved it by 
exact numerical diagonalisation of Hamiltonian Q rep- 
resented in real-space. We consider a problem of N sites 
(up to N = 1728) with one special site, and we calculate 
all the eigenvalues and the total energy. A key question 
is whether bound states with energies lower than that al- 
lowed by band-structure can appear near the defect. We 
find that they do in several cases, and when their bind- 
ing energy exceeds the gap, it signals an instability of the 
periodic phase. 

A. Instability of the Ferromagnetic Insulating 
Phase when g/t is decreased 

We consider first the FI-CO phase at half-doping (pic- 
tured in Fig. |2J). It is stable for very small Jaf and large 
g/t (see the phase diagram in Fig. Out of the two 
sites in the unit-cell, one site is distorted with a distortion 
orientation that favours the d x 2_ y 2 orbital. 

We now consider the problem when one introduces a 
single site JT defect: the amplitude of the distortion Q of 
the FI-CO phase is maintained at N/2 — 1 sites except at 
one site where the distortion is now reduced to Q — Qd- 
Qd takes all values from (no defect) to Q (the lattice 
distortion has been completely removed on this site). The 
excess energy of such a state is given by: 

E-E = E el (N 0: Q, Q d ) - E el (N 0: Q,0) 



+ \k{Q - Q d ? - \KQ 2 (6) 

Here E and E el (N ,Q,Q) [N = N/2 is the number of 
electrons] are the total and electronic ground state ener- 
gies of the optimal FI distorted phase, obtained as a func- 
tion of g/t by minimising with respect to Q as discussed 
in the previous section. E e i(No,Q,Qd) is the electronic 
ground state energy of the defective state. One expects 
a gain in lattice energy and a loss in electronic energy, 
because one energy level has been raised at one site. To 
evaluate the latter, we first solve numerically the problem 
of the one electron eigenvalues in the presence of the ex- 
tra single site potential for a finite-size system. Then we 
calculate E e i(N ,Q,Qd) by filling the N/2 lowest one- 
electron levels. We have so far considered 3d systems 
with up to N — 1728 sites. 

Fig. |S] shows the energy E — E plotted vs. Qd 
for different values of g/t. We have checked that finite 
size effects are negligible (the curves corresponding to 
N = 216, 1000, 1728 are given for g/t = 6.7 in Fig. 0). 
For large g/t, the energy is positive but there is a lo- 
cal minimum at large Qd which can be described as a 
particle-hole excitation with reduced distortion on one 
site. When g/t decreases, this excitation softens and 
vanishes at g c /t ~ 6.8. We believe that this signals the 
onset of a new phase where such defects are energetically 
favourable and thus proliferate in the system. 
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FIG. 8: Energy change when a single JT defect is introduced 
in the FI-CO phase, vs. Qd- Q — Qd is the JT distortion on a 
defect site; all the other occupied sites having the same distor- 
tion Q. We see that there is an excitation with Qd ~ Q that 
softens when g/t decreases. The excitation corresponds to a 
band particle-hole excitation with the removal of a lattice dis- 
tortion of one site, while the Qd = minimum is the polaron. 
The softening for g c /t ~ 6.8 signals a phase transition with 
proliferation of mobile electrons and defects. Finite-size ef- 
fects are small and shown for g/t = 6.7 (TV = 216, 1000, 1728). 

From the calculation of the energy levels in the pres- 
ence of the defect, we find that there is no bound state 
within the gap for the Qd that minimises the energy. The 
electron occupies a higher-energy band-like state and is 
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mobile. The instability therefore corresponds to the en- 
ergy of this mobile electron crossing the chemical poten- 
tial (i.e., —Ejt)- This suggests that the proliferation of 
the defects leads to the conversion of some small fraction 
of the localised electrons into mobile electrons moving 
on weakly distorted sites, resulting in a metallic phase. 
Such a state would not be accessible in the minimisation 
procedure of section ITU ( which has a maximal unit cell of 
8 sites) even if the defect sites were to arrange themselves 
in a super-lattice. 

The situation is rather similar to that described by 
Ramakrishnan et al. at the metal-insulator transition in 
hole doped manganites in the orbital-liquid regime^ ex- 
cept that we have here a calculation in the context of a 
microscopic model that explicitly suggests such a picture 
even when g/t is not very large. Thus we can identify 
the high-energy mobile electrons as the broad-band Mike 
electrons of ref . , and the low-energy localised states 
as the I polarons. In the present context, at half-doping 
and above the transition (<? > g c ), all the electrons oc- 
cupy the I states, which form a regular checker-board 
array (Fig. EJ). The sites are singly occupied and U 
does not play a crucial role. This is no longer the case 
below the transition when we start to transfer some elec- 
trons from the £ states to the b states. The b states are 
delocalised over the empty sites but also visit the sites 
occupied with £ electrons. Double occupancies become 
inevitable and U has to be taken into account in order to 
determine accurately the properties of the metallic state. 
The question of what kind of new metallic state arises 
for JT couplings just below the instability is clearly in- 
teresting. The mobile electrons, for instance, may be 
able to destroy the orbital and charge order. While a 
study of such issues is beyond the scope of the present 
article, the above results suggest than one should add 
mobile electrons to the strong-coupling Hamiltonian JSJ, 
in order to describe metallic phases with possible partial 
orbital/charge order: 

H = —Ejt fin + J^Si.Sj 

i <ij> 

~ E ( S « ■ S * + Ml - n*)C& + (* ~ i)] 

<i,j> 

- E C^+E^W* ( ? ) 

<i,j> i 

where the orbital index a of the mobile b electrons takes 
both values on the undistorted sites, but is constrained 
to be orthogonal to the I orbital on the occupied sites, 
and the other quantities have the same meanings as in 
eq. 0. For infinite U the mobile electrons can- not hop 
to the I sites at all, and the last pair of terms can sim- 
ply be replaced by - T,<i,j> CALAM 1 ~ n a){ 1 ~ n tj)- 
This Hamiltonian needs to be studied in a framework 
that can handle the strong interaction effects, such as 
the dynamical mean-field theory, in a similar way as was 
done before for the orbital-liquid stated Note further- 



more that the above Hamiltonian does not include i — b 
hybridisation effects, which must be included in order to 
describe properties sensitive to I — b coherence which can 
develop at sufficiently low temperatures in the metallic 
phases^ It is straightforward to generalise the Hamilto- 
nian to include these effects, as well as cooperative JT 
effects. 



B. Instability of the CE Phase upon Doping 

In the band picture of the CE phas o 28 i 52 , doping with 
electrons, corresponding to x < 1/2 (resp. holes, corre- 
sponding to x > 1/2) provides mobile carriers in the con- 
duction (resp. valence) band. In either case, the system 
will be metallic. This is contrary to experiment, where, 
in most cases, the system remains insulating for x > 1/2, 
but typically becomes metallic quickly for x < 1/2. The 
transition to the ferromagnetic metal for x < 1/2 has 
been described as being due to the crossing of the energies 
of the CE and ferromagnetic metallic states^ The tran- 
sition is then naturally first-order. However, as discussed 
in ref. |33| even for x < 1/2 a simple band picture of the 
ferromagnetic metallic state is severely limited. Apart 
from that, the band picture fails to describe the insulat- 
ing character of the regime x > 1/2 and the particle-hole 
asymmetry around x — 1/2. We discuss this issue next. 

It was pointed out a long time ago by de Gennes^, 
in the context of slightly doped LaMn03, that adding 
carriers to the antiferromagnetic phase of LaMn03 may 
favour canted structures. The qualitative argument is 
that at small concentration the carriers gain kinetic en- 
ergy which is linear in the canting angle whereas the loss 
of magnetic energy is quadratic in the canting angle. By 
the same token, adding carriers to the CE phase should 
lead to canting of the core-spins. As such phases inter- 
polate between the CE and FM phases, the transition to 
ferromagnetism should be naively second-order. 

In view of this, we have calculated the energy of homo- 
geneous CE canted phases (defined in Fig. EJ) for differ- 
ent carrier concentrations on either side of x = 1/2 (i.e., 
retaining the 8-sublattice periodic structure even when 
x 1/2). We find that canting is favourable for adding 
electrons to the half-doped system but not for adding 
holes, as de Gennes's general argument is valid only for 
very small carrier concentration and breaks down quickly 
on the hole side, due to special features of the CE state. 
We have calculated, in addition, the energy in the pres- 
ence of a single-site defect in the JT distortion as in the 
previous subsection. We find that, when g/t is sizeable, 
canting is in competition with self-trapping of the carri- 
ers in JT defects. We find that the energy gain due to 
trapping is linear in the carrier concentration (and thus 
dominates at low concentration) whereas it is quadratic 
for the canting. For intermediate values of g/t, this re- 
sults in a first-order transition to a canted metallic phase 
when electrons are added to the half-doped system (i.e., 
for x < 1/2), and the persistence of a CE-type phase 
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with self-trapped carriers when holes are added (i.e., for 
x > 1/2). 



1. Canted Phases 

We first consider the total energy of a homogeneous 
canted CE phase as a function of the canting angle (j> (see 
Fig. [5] for definition) for different values of doping (Fig. 
1101 top). We take the JT distortions that optimise the 
energy (as discussed in section|nJ, except that we neglect 
the small distortions on the corner sites for simplicity. 53 
When adding electrons (x < 0.5), the energy of the sys- 
tem is lowered by a finite canting angle (Fig. top). 
The higher the concentration the higher the canting an- 
gle. When adding holes (x > 0.5), however, it turns out 
that the system prefers = 0. 




FIG. 9: (color online). Canted CE phase with canting angle 
cj> (CE-M-C). The electrons can now hop from chain to chain. 
The original Id zig-zag chains are marked with thicker lines. 

We can understand these results by considering in more 
detail the band structure of the CE phasei 28 i 52 ' 54 As dis- 
cussed earlier for x = 1/2 (section^), the bridge sites are 



orbitally ordered in the (3x 



! )/(3y 2 



! ) pattern. 28 



If the Jahn- Teller distortions occur in such a way as to 



favour precisely the alternation of (3x 



! )/(3y 2 



orbitals, the energy of the system is further lowered. The 
band structure has four dispersive bands and four non- 
dispersive bands (two at zero a nd two at finite energy), 
as shown in Fig. ^2 (see Refs. 28 52 54])and described 
by: 



e 3,4 



= 0, 



- E j T ± \J Ej T +t 2 (2 + cos q a ) 
-E JT ± yfE 2 JT +P(2-cosq a ) 

£5,6 = 2Ejt 



(8) 

(9) 
(10) 



where t = At/3 and Ejt = gQ/2, and the wave-vector 
q a , which takes values in the reduced Brillouin zone 
[— tt/2, 7r/2], is parallel to the chain direction. The band 
structure is analogous to that of Rcf. with a 

charge ordering coming from the Jahn- Teller distortions, 
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FIG. 10: Top. Energies of the canted phases as function of 
the canting angle, for different concentrations (g/t = 5 and 
Q = Qopt = 0.35). For x < 0.5, there is a finite angle that 
minimises the energy. For x > 0.5, the angle is zero and the 
CE state is stable. Bottom: Comparison of the energies of 
the different phases. The canting angle is chosen such as to 
minimise the energy as in the top figure. 



as in Ref. |5J|, but determined self-consistently. At 
x = 1/2, the lowest band is completely filled. 

The zero energy band is made up entirely of the states 
from the corner sites (see appendix A). The charge gap, 
from the top of the filled valence bands to the zero energy 
states, is given by: 



E 



JT 



Ej T +t* 



(11) 



When the core spins are canted away from = 0, the 
degeneracy of the zero energy states of the zig-zag chains 
is lifted, and they form bands which disperse. To first 
order in the canting angle, the dispersion arises from the 



coupling of the (3z 



orbitals at the corner sites along 



the z-direction, and is given by: 



e°(q a ,q z ) = -At<p 



1 + COS q a 

2 + cos q a 



cosq z 



(12) 
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FIG. 11: Band structure of the Id zig-zag chains (solid lines) 
with Jahn- Teller distortions (Ejr/t — 0.875), given by eqs. 
©, © and itTUli 28 ! 52 ! 54 At x = 1/2, only the lowest band 
is filled. Also shown is the splitting of the zero-energy states 
(dashed lines) when the core spins are tilted away by an angle 
<j> (whence their dispersion becomes 2-dimensional) . 

Added electrons (with respect to the reference state 
x = 1/2) occupy the bottom of this new band (prop- 
erly folded in the reduced Brillouin zone, Fig. Illfl at 
q, q z ~ 0. Each electron therefore gains an energy —f3 e <fit 
with (3 e = 8/3, which is completely independent of g. On 
the other hand, there is an energy loss K(j> 2 per site, where 
k is the effective spin stiffness towards canting. There are 
two contributions to k. One comes from the direct super- 
exchange. In the CE phase with cf> = 0, two of the neigh- 
bour spins are parallel and the other four neighbour spins 
are antiparallel. This gives a contribution AJafS 2 4> 2 to 
the energy. There is also a contribution from the double 
exchange, which we calculate numerically by calculating 
the change in the total kinetic energy as function of <f>. 
From this we extract a stiffness in the limit of small </>, 
~~ K e = f(g/t, 5c) where 5c = 1/2 — x is the filling fraction 
in excess of that at half-doping, so that k = AJafS 2 — K e . 

The canting angle is then given by minimising the ex- 
cess energy per site, E — E — n<j) 2 — f} e t5c<j>, and is given 
by 4> — (3 e t5c/2n. Note that the linear dependence in 
5c is valid only for sufficiently small 5c whence the ad- 
ditional electrons occupy states near the bottom of the 
band. There are in fact N/8 states with energy smaller 
than and N/8 states with energy greater. So the linear 
dependence is expected to be reasonable for 5c <C 1/8. 
The total minimised excess energy is then 



E-E = -^(5c) 2 : (13) 

quadratic in the electron concentration (Eq is the energy 
of the CE phase at x = 1/2). The canted CE phase thus 
has a lower energy than that of the CE phase by this 
amount, whenever 5c > 0. This explains the quadratic 
behaviour in 5c found numerically in the previous para- 



graph (Fig. ^3 bottom). When 5c increases further, the 
canting angle increases until eventually the system be- 
comes fully ferromagnetic via a second-order phase tran- 
sition. Within the above picture, as soon as 5c > 0, 
the system is metallic because the additional electrons 
occupy the dispersing conduction band. But note that, 
since the dispersion is mainly along the z and the chain 
directions to leading order in the canting angle, the met- 
alicity generated by such a mechanism would be highly 
anisotropic. 

Surprisingly, the situation turns out very differently 
when excess holes are added to the half-doped system. 
Naively, one might have thought that canting the spins 
will push up the states at the top of the lower band and 
that the added holes will thus gain energy, similarly to 
the case of electrons. But according to the numerical 
results of Fig. ^| this is not what happens. The energy 
is minimal at (f> = 0, i.e. when the spins are not canted. 
The reason is that the energy gain for the whole system 
is not linear in the number of additional holes, except for 
extremely small 5c. This is even more emphatically and 
dramatically evident from Fig. which shows that the 
energy gained by canting for a small fixed canting angle 
(cf) = 0.05) as a function of doping, is extremely particle- 
hole asymmetric. While the gain is indeed linear in the 
deviation from half doping on either side for the smallest 
values of \5c\, it drops quickly and substantially below 
this, and stays small for x > 1/2, i.e., for added holes. In 
contrast, when electrons are added, the energy gain from 
canting remains large for large \5x\. (We emphasise that 
the above results are only for the canted CE state for any 
x; hence they are most meaningful near x = 1/2.) 

The reason for this asymmetry is the difference in the 
dispersion of the electron and hole bands in the canted 
CE phase. As discussed in detail above, the conduction 
band in the canted CE phase (into which electrons get 
added) disperses only in two directions, and hence the 
corresponding density of states is constant, over a band- 
width of order t<j>. Hence the energy gain from canting 
remains substantial up to half occupancy of the conduc- 
tion band, corresponding to 5c = 0.25, beyond which the 
energy gain starts reducing because of cancellations, as 
is clear from Fig. ^] (dashed line). In contrast, the va- 
lence band in the canted CE phase which accommodates 
the holes (and which we have not discussed in detail), 
disperses in all three directions, resulting in a density 
of states which starts from zero at the top of the band. 
Typically when the change in the chemical potential is 
smaller than tip, the additional holes occupy the higher 
energy states of the new band and the energy gain is lin- 
ear in the number of additional carriers. On the other 
hand, when the chemical potential is of order tip, both 
sides of the new band are occupied and there is no en- 
ergy gain. The condition is expressed as 5c < t(pp(tp) 
where p(ef) is the density of states at the Fermi energy 
for the system with the concentration c + 5c. When the 
density of states is small, the range within which the en- 
ergy gain is linear in the number of carriers is small. In 
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this range and out of it, the canting angle is then very 
small because it does not scale with the number of carri- 
ers anymore, as is confirmed by Fig. El 
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FIG. 12: Energy gained by canting of the CE phase for a fixed 
canting angle of <j> — 0.05 as a function of doping away from 
half doping, showing the particle-hole asymmetry. 



2. Trapping of Added Carriers in the CE Phase 

Trapping of Added Electrons: We next explore the 
competition between the above process and the possible 
trapping of the added electrons in lattice distortions by 
allowing for an additional lattice distortion on one spe- 
cial site and iVo + 1 electrons (No = N/2). In other words 
we look at the cost or gain of energy involved in trapping 
one additional electron in the CE state. Since the added 
electron goes into a band made up from states belong- 
ing to the corner sites which is originally undistorted, 
we choose the special site with the added distortion to 
be one of the corner sites, and of such an orientation as 
to lower the energy of the (3z 2 — r 2 ) orbital (Fig. I13fl . 
As in llll Al translation invariance is now broken, and we 
find the one-electron energy eigenvalues by diagonalizing 
exactly the problem on large lattices, as a function of 
the strength of the additional distortion, Qd- We show 
the energy levels in Fig. ^0 left. In addition to the band 
states (black areas), we also find a couple of bound-states 
within the gap (Fig. 1131 left). 

The total excess energy including the cost of elastic 
energy to create such a defect in the lattice is calculated 
by filling up the energy levels with No + 1 electrons: 



E - E a = E el (N + l,Q d )-E el (N + 1,0) + -KQ 2 d (14) 

The energy can be viewed as the sum of three different 
contributions, each of which is separately shown in Fig. 
IT31 there is the electronic energy gain for the electron 
bound to the defect, the scattering energy for the elec- 
trons that are scattered by the defect (both of which are 
contained in E e i(No + 1, Qd)) and the elastic energy cost 
for creating such a distortion. When all are put together, 
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FIG. 13: (color online). Definition of the defect (top), that 
is an additional distortion of strength Qd on a single corner 
site ((3z 2 — r 2 ) orbital). Energy levels vs. Qd (bottom left). 
At Qd = 0, there is no defect and the band structure is that 
of Fig. 1111 At Qd > 0, bound states appear within the 
gap. Total excess energy of the system with No + 1 electrons 
(bottom right). This can be seen as the sum of the energy 
of the bound state Eb s , the elastic energy E e i ast and defines 
the scattering energy E scat t- g/t = 5, Q = Q op t = 0.35 
(Ejt = 0.875). 



it turns out that it is favourable to create a defect with 
an energy gain, E JT . For instance, in Fig. 1131 when 
g/t = 5, the energy gain is E JT = 0.75t ~ 0.86E JT , 
slightly smaller than Ejt due to band structure effects. 
The problem can not, indeed, be reduced to polaron for- 
mation purely with the states of the conduction band by 
ignoring the valence band (that would give E JT = Ejt)- 
The scattering of electrons by the defect and the level 
repulsion between the two bands (Fig. I13fl accounts for 
the reduced value of E e JT found above. 

For a small concentration of additional electrons, 5c, 
the energy gain with respect to the CE phase, given by 



E-Eo 



-Ej T Sc, 



(15) 



is linear with the concentration of additional electrons. 
Here Eq is the energy of the perfect CE phase with addi- 
tional electrons occupying the (undistorted) zero-energy 
band states. 

When we compare eq. (TPol) and eq. (JT^J, it is clear 
that the linear dependence (fPo)) gains over the quadratic 
dependence (|13f) of the canted phases at small doping. At 
small doping then, the system is insulating and the addi- 
tional carriers are trapped in inhomogeneous lattice dis- 
tortions. The resulting phase is interesting in that it pos- 
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sesses some self-induced disorder. Random self-trapping 
of additional carriers has also been reported in the adi- 
abatic spinless Holstein model, with a concentration of 
carriers close to one electron for two sites ^ that is simi- 
lar to the present case. By including the double exchange, 
we have shown that such doping-induced inhomogeneous 
states are in fact quite general and more stable than 
the canted states because of the linear energy gain we 
have found at small concentration. It has some similari- 
ties with the inhomogeneous states with metallic droplets 
found in a simple model with charge-ordering near half- 
filling^ We note that we have considered here the sim- 
plest polaronic state with fully localised electrons on the 
defect sites. It would be interesting to study whether 
more complex defects, involving for instance distortions 
of the magnetic structure on the neighbouring sites (mag- 
netic polaron) could arise near x ~ 1/2, as suggested for 
x ~ (ref. 0) or x ~ 1 (ref. [13). 

Trapping of Added Holes: Similarly, it is also 
favourable to trap added holes in lattice distortions. We 
consider the analogue of the previous problem with one- 
less distortion and N — 1 electrons. The removal of the 
distortion on one site creates again a defect. The energy 
levels and the excess energy calculated numerically by 
filling the energy levels with N n — 1 electrons, 



E-E = E el (N - 1, Q, Q d ) - E el (N - 1, Q,0) 

+\K{Q-Q d f -\KQ 2 (16) 

are given in Fig. [21 We find that it is favourable to trap 
the additional hole onto the bound state that appears 
within the gap just above the lowest band (Fig. [21 left) 
as soon as g/t > 4. The energy gain, Ej T is smaller than 
E JT and is 0.31t ~ 0.35.E 7T for g/t = 5 (Fig. EJ). For 
the whole system the energy gain is —Ej T 5c. It is then 
favourable to trap the holes and the system is insulating. 



3. Competition between the two phases and Phase Diagram 

By comparing the energies of the canted state and the 
defective state, (for instance, as shown for g/t = 5 in 
Fig llUfl we can arrive at a phase diagram in g/t~x plane 
near half doping. As discussed above, when electrons are 
added to the CE phase, they are trapped on corner sites 
with newly generated Jahn- Teller distortions if their con- 
centration is within [1/2 — 6c cr it, 1/2] fFig. I15|) . These 
JT distortions are oriented in the ^-direction, so as to 
favour the occupancy of the (3z 2 — r 2 ) orbital. The mag- 
netic structure remains of CE type. For x = 1/2 — 5c cr it 
there is a transition to a metallic state with canted spins. 
In the present approach there is a finite canting angle at 
the transition, so that the transition is first-order. When 
the concentration increases further the canting angle be- 
comes larger and larger. At small g/t we find that there 
is a first-order transition line between a canted state with 
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FIG. 14: (color online). Definition of the defect on a single 
bridge site with distortion decreased to Q — Qd (Qd = cor- 
responds to a fully distorted site as in the original structure) . 
Energy levels vs. Qd (bottom left) and energies for No — 1 
electrons (bottom right), g/t — 5. 



small canting angle and a highly canted state that ends 
by a critical point (Fig. I15fl . The system becomes even- 
tually fully ferromagnetic. 

On the other hand, when holes are added to the CE 
phase, the holes get trapped on bridge sites with lat- 
tice distortions removed when g/t > 4 (and the system 
remains insulating, CE-T) or move freely in the lowest 
band for g/t ;$ 4 (and the system is metallic, CE-M) 
[Fig. I15| . There is no competition with canted states 
in this case because, as discussed in IIII B ll the density 
of states near the top of the valence band is not large 
enough to provide sizeable energy gains. 

We emphasise that this phase diagram is based on an 
instability calculation, and can be expected to be accu- 
rate only when x is sufficiently close to x = 1/2. Away 
from it, the defects start to interact and other phases 

12 20 21 

may appear iia*aiaai 

We note that in all our discussions so far we have ne- 
glected the disorder that is present in the real systems. 
In fact, even at half-doping the arrangement of the A site 
ions (La or Ca) is disordered. This disorder causes the lo- 
calisation of the one electron states (near the band-edges 
in 3-d) leading to increased stabilisation of the insulat- 
ing properties near half-doping, and may also have con- 
sequences on the local magnetic structure as emphasised 
by the idea of ferronsi^ 
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FIG. 15: Phase diagram as a function of the model parameter 
g/t for hole doping concentration x close to 1/2 (JafS 2 /t = 
0.15). CE-M-C denotes a metallic and canted CE phase. CE- 
T refers to the CE phase with added carriers self-trapped in 
JT defects^ CE-M is the CE phase with added holes that is 
metallic (no canting nor trapping of the holes). C is a critical 
point ending a first-order transition line between two CE-M-C 
states with differing canting angles that arise for small g/t. 
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FIG. 16: First-order transition of the CE phase to a canted 
state when an external magnetic-field is applied to the system 
at g = 0. The threshold field that induces the first-order 
transition is given by g/i_g_ff c = 0.0KM. The system is metallic 
beyond H c . 
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C. Instability of the CE Phase in a Magnetic-Field 

Finally, in this subsection we discuss the instabilities 
of the CE phase in a magnetic-field due to the combined 
effects of canting of spins towards the direction of the 
field by a canting angle (j> (as pictured in Fig. [5} and 
modification of the JT distortions. 

Since a band opens out of the zero energy states upon 
canting (dashed lines in Fig. Illf) . the system will "melt" 
into a metallic phase when the corresponding charge gap 
closes. Since this can happen for a canting angle less 
than tt/2, the transition will typically not be to the fully 
polarised ferromagnetic state. 

In Fig. 1161 we show the energy versus the canting angle 
for different fields at g/t = 0. There is a small optimal 
canting angle that minimises the energy for small fields, 
but the system continues to be insulating. When the field 
increases further there is a first-order transition to a state 
with a finite (and large) canting angle, which is metal- 
lic. This is reminiscent of the first-order transition line 
between two canted states when additional electrons are 
added (section flll Bl and Fig. I15|l . We note that the first- 
order transition field is very small, gfisH c = O.OlOi, not 
only because, with JafS 2 = 0.15i the system is close to 
the phase boundary, but also because the canted phases 
reduce the critical field considerably. Previous works that 
have assumed that the insulator-metal transitions involve 
the fully ferromagnetic state2£i2L2£ would predict a tran- 
sition field g/is_ff c = Q.lAt for the same Jaf we have used 
above. We discuss in detail next as to how these features 
change with respect to turning on g (and keeping the 



To start with we take the CE phase with the optimal 
lattice distortions found in section [D] for different values 
of g/t. First, we freeze the distortions for all magnetic- 
fields and find the optimal canting angle that minimises 
the energy as a function of the magnetic-field, for non- 
zero values of g/t (Fig. I17|) . The jumps in the canting 
angle correspond to first-order transitions between states 
with small and large canting angles. When g/t is in- 
creased, the transition fields and the fields at which the 
fully ferromagnetic state is reached shift rather quickly to 
larger values. The transition field becomes larger simply 
because, with increasing g/t (and for the same Jaf), the 
system moves further and further away from the phase 
boundary with the ferromagnetic phase (see the phase 
boundary in Fig. [2J) . 

We next show that the threshold fields for the "melt- 
ing" of the CE phase are reduced by taking into account 
the effect of the magnetic-field on the lattice distortions 
themselves. One should expect this, since in section ITT1 it 
was shown that the ferromagnetic phases remains undis- 
torted up to g/t — 5. It is made explicit in Fig. 1181 
where the energy of the undistorted highly-canted state 
(solid line), which merge into the fully ferromagnetic 
state (dashed line) at higher fields, crosses that of the 
distorted canted CE states. These transitions occur at 
smaller transition fields (see Fig. I18|) compared with the 
transition fields we considered in the previous paragraph 
(the latter correspond to the cusps, visible in the curves 
g/t = 2,3 in Fig. I18|l . The new magnetisation curves 
are shown in Fig. ^] top. We then calculate the band 
structure as a function of the magnetic field and extract 



16 





1 1 1 1 




- y 






/ 






/y^ g/t=o-5 - 




J AF S - 
1 1 1 1 


0.15t 



0.1 0.2 0.3 0.4 0.5 
g^ B H/t 

FIG. 17: Optimal canting angle vs. external magnetic-field 
for different values of g/t (from up to 5) keeping the same 
distortions as at H — (CE phase) for all magnetic-fields. 
There are first-order transitions between canted states. As 
we will show, there are other transitions, involving modified 
JT distortions, that preempt those shown in this figure (see 
below, Fig. [T51 . 



the gap (Fig. 1191 bottom). The jumps in the magneti- 
sation and in the gaps turn out to be simultaneous, thus 
indicating that the transitions correspond to insulator- 
metal transitions. The transition fields g/is-ff c vary in 
the range O.Oli - 0.2t for J A fS 2 = 0.15i. The smaller 
the g/t the smaller the transition field at fixed Jaf- 
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FIG. 18: Crossing of the energies of the distorted canted 
CE phases (solid lines, with g/t = 2,3,4,5) with that of 
the undistorted highly-canted phase (solid line corresponding 
to g/t = 0), or the fully ferromagnetic phase (dashed line) 
to which it merges for large fields. The cusps visible in the 
g/t = 2,3 curves at higher fields correspond to the first-order 
transitions described in Fig. 1171 But these are preempted by 
the transitions to the undistorted phase. 

The first-order transitions discussed above were ob- 
tained by considering the crossing of two solutions, 
namely the distorted canted CE phase (with the distor- 



tions frozen at their H = values) and the undistorted 
FM phase. It is possible, in principle, that intermediate 
phases with intermediate distortions appear. To rule out 
this possibility, we have performed the full optimisation 
in presence of the external magnetic-field for g/t = 4,5 
and found that the decrease in the distortions is less than 
6% up to the transition to the undistorted ferromagnetic 
phase. This validates our approximation of using frozen 
distortions in the canted CE phase up to the first-order 
transition. 




gH B H/t 

FIG. 19: Magnetisation or canting angle (top) and charge 
gap (bottom) vs. external magnetic-field for different values 
of g/t = 1 - 5 (J A FS 2 /t = 0.15). There are clear first-order 
transitions to an undistorted canted metallic or fully ferro- 
magnetic metallic phase. 

We can roughly understand some of the above results 
in terms of an expansion of the physical quantities in 
the canting angle, before the first-order transition takes 
place. For simplicity we neglect the distortions on the 
corner sites of the CE state, so that the zero-energy states 
are degenerate to start with. The degeneracy is then 
lifted by the canting and the reduced charge gap is given 
at first order in the canting angle by: 

A(0) = A - f3 e t<t> (17) 

where A is the gap of the uncanted CE phase, given 
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by equation (|llf> . and (3 e = 8/3 a numerical coefficient 
(from eq. 11211 . The canting angle itself is determined by 
minimising the magnetic energy per spin due to canting, 
given for small 4> by: 

E = ncj} 2 - gfisSHcj) (18) 

where k = 4JafS 2 — n e is the spin stiffness we have re- 
ferred to before fsection llll Bp . The energy is minimal for 
4>opt{H) = gfiBSH/2n, which describes the linear regime 
of the curves in Fig. El before the first-order transitions 
to the highly canted states take place. The corresponding 
charge gap, given by 

A mm = A(cj> opt ) = A - (^j gn B SH, (19) 

also decreases linearly with field, but is not a very good 
approximation to the results of Fig. El bottom. The 
difference comes from the small distortions on the corner 
sites which we neglected in deriving (|17J) by assuming 
degenerate zero-energy states. With small distortions on 
the corner sites, the zero-energy states are no longer de- 
generate, and that would modify the expression for the 
gap J32|). 

In any case, one can give an upper bound for the criti- 
cal magnetic-field by extrapolating the linear expression 
(|TT)j) to zero. This is an upper bound because the actual 
transition occurs before the gap get completely closed as 
is clear from Fig. El We have 

gHBSH c = ^A = (8 Jaf S 2 - 2n e ) A (20) 

Pet Pet 

It is clear from this expression that in order to close a gap 
A of order Ejt or t, we do not need a magnetic-field of or- 
der A thanks to the reducing factor njt. n/t is small, first 
because the spins are easy to polarise on energy scales 
that have nothing to do with the charge scales, as in 
the standard field-induced insulator-metal transition of a 
spin-density wave. It is interesting to note that the large 
field gfiBSH c o = 8JafS 2 is the critical field to align anti- 
ferromagnetic spins coupled only by Jaf- The real tran- 
sition fields are substantially reduced compared to this by 
the double exchange included in K e and the factor A/ [3 e t 
which describes how fast the gap closes with the cant- 
ing angle. For example, the strength of the transitions 
fields of Fig. El are i n the range g^sH c ~ O.li — 0.25i. 
With t ~ 0.2 eV, we have gfJ, B H c - 20 - 50 meV, or 
H c ~ 140 — 350 Tesla, which are still much too large. 
We discuss this discrepancy below. We now consider the 
other situation with g/t > 5.0. 

2. Transition to an Inhomogeneous State 

For g/t > 6.8 we have seen that even the ferromag- 
netic phases are insulating. Hence it is to be expected 
that no insulator-metal transition can take place in this 



regime. Instead, the CE phase makes a transition to 
the ferromagnetic insulating phase for sufficiently large 
fields, as is clear from the large g/t limit phase diagram 
discussed in section lll Dl For g/t < 6.8, however, we have 
shown earlier that it is favourable to create defects and 
mobile electrons out of the FI-CO phase (section llll All . 
We expect therefore that similar kinds of phases will be 
favoured by a magnetic-field as well. We explore this 
issue further below, taking g/t = 6 as an example. 

We first compare the energy of the canted distorted 
CE phase [with the usual (3x 2 — r 2 )/(3y 2 — r 2 ) orbital 
ordering] with that of a canted phase with alternate sites 
distorted in such a way as to favour the (x 2 — y 2 ) orbital 
ordering. The latter is indeed a lower energy state for 
large g/t when the spins are fully aligned, compared with 
the undistorted phase. We see in the inset of Fig. I20lthat 
the two energies cross for g^sH/t ~ 0.43. As shown in 
Fig. 1201 at this field the canting angle that minimises 
the energy jumps from ~ 0.8 to 7r/2 corresponding to 
a transition to the fully ferromagnetic state (which is 
metallic at g/t = 6.0, see Fig. [SJl with (x 2 — y 2 ) orbital 
ordering. 

We now consider the possibility of creating a defect 
in the lattice distortion pattern, similarly to what we 
discussed in section 1111 Al which might lead to a lower 
energy according to what we found there. We start with 
the [x 2 — y 2 ) phase with canting angle cf>, and reduce the 
distortion of one of the distorted sites to Q — Qd instead 
of Q, so that when Qd = Q, the distortion is completely 
removed. In Fig. [21 we show the total energy E — Eq 
(where Eq is the energy of the homogeneous phase, cor- 
responding to Qd = 0) as function of Qd for different 
canting angles near the transition. Clearly, there is a 
minimum at Qd — Q {Q = 0.39 at g/t = 6.0) that corre- 
sponds to an excitation where we remove a distortion and 
create a particle-hole excitation. This minimum becomes 
soft at <f> ~ 0.9, thus signalling the onset of a transition to 
a phase where such defects are energetically favourable. 
Note that this phase can not correspond to the undis- 
torted phase (i.e., to removing the distortions on all the 
sites) because the undistorted canted phase is higher in 
energy as is clear from Fig. 1201 fshort-dashed lines). 

To sum up the discussion above, when a magnetic-field 
is switched on, the core spins cant towards the direction 
of the field. A band opens out of the zero-energy states of 
the corner sites of the CE phase and the gap gets reduced. 
At a threshold value of the magnetic-field, there is an in- 
stability toward transferring electrons from the localised 
states (or lower energy band) to the mobile states. For 
field values close to this threshold, further energy gain is 
possible when some Jahn- Teller distortions are removed 
and defects are created. When the density increases, 
the missing distortions start to play a more important 
role and the one-defect approach we have developed here 
breaks down. This scenario for the magnetic-field in- 
duced metallic phase is very close to the description of 
the metallic phase in the (x < 0.5) regime as discussed 
in Ref. |33fl . 
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FIG. 20: Magnetisation or canting angle (top) vs. external 
magnetic-field for g/t = 6 (JafS 2 /t = 0.15). In the inset the 
energies of the undistorted canted phase (dotted line), the 
canted phase with distortions so as to favour (x 2 — y 2 ) orbital 
ordering (dashed line) and the CE canted phase (solid line) are 
given. The latter two cross for gfisH/t around 0.43, where 
the first-order transition takes place. Above this field the 
phase is fully ferromagnetic. Bottom: The resulting Charge 
Gap is shown as a function of qhbH /t 



Clearly, for all the values of g/t we have considered, 
there are rich transitions from the CE phase to the fer- 
romagnetic phases as a function of the magnetic-field. 
The lattice distortions play a crucial role in converting a 
transition that would be naturally second-order because 
of the progressive canting of the spins into a first-order 
transition. 

The strength of the magnetic-field at which the tran- 
sition occurs, although much reduced compared to naive 
estimates as discussed above, is still much too large com- 
pared to experiment, especially if we consider large val- 
ues of g/t. However, the large distortions observed ex- 
perimentally are consistent with relatively large values 
of g/t. How does one reconcile these two? As we have 
seen, a crucial ingredient that determines the strength of 
the transition field is the charge gap. We can reasonably 
argue that the gap is overestimated in the present ap- 




FIG. 21: Energy for creating a single-site defect in the JT 
distortions with amplitude Qd for various canting angles in 
the FI phase. Such an excitation becomes soft at <j> — 0.95 
and the defects proliferate. The transition corresponds to an 
instability towards a metallic phase with defects. 



proach. On the experimental side, the gap is not of the 
order of t but much smaller, typically 5 times smaller for 
Ndx ^Ca^MnC^^ This is also clear from the tempera- 
ture at which the charge-order transition takes place. If 
we use the estimates of the gap from experiments in the 
expression l|2U[) . the transition field comes out to be much 
smaller, of order 30 T, much closer to the transition fields 
of that compound. 

It is clearly important, therefore, to improve the the- 
ory presented above to generate a more accurate esti- 
mate for the charge gap in the CE phase. For this we 
need to include at least three sources of correction to the 
gap estimate above. Namely the finiteness of the Hund's 
coupling Jh, the cooperative nature of the JT distor- 
tions, and small second neighbour hoppings. When Jh 
is finite, the electrons can hop even between sites which 
have anti-aligned core spins. This opens a band out of the 
zero energy states and reduces the gap. Second, because 
distortions of neighbouring sites are coupled, a distor- 
tion on a site imposes a distortion on the neighbouring 
site that lowers the energy level, thereby reducing the 
gap. Finally, second neighbour hopping allows the corner 
site electrons to become mobile (even in the case of infi- 
nite Jh), and overall increases the bandwidths of all the 
bands, hence reducing the charge gap (see, for example, 
ref. [54()- We believe that, while the theory presented 
above clarifies how the smallness of the transition field 
arises, a more elaborate theory including these three ef- 
fects which reduce the gap is required to obtain a precise 
estimate of the transition field. 



IV. CONCLUDING REMARKS 

In conclusion, we have first confirmed various periodic 
phases of the phase diagram of half-doped manganites 
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(Fig. 0) , by optimising the lattice and magnetic energies 
in the thermodynamic limit. We have thus provided ex- 
plicit calculations for the Jahn- Teller distortions, charge 
and orbital order-parameters in the various stable phases 
by exploring the phase space in an unbiased way, albeit 
limited by the 8-site unit-cell. 

It is interesting to discuss, in the context of the phase 
diagram of Fig. the strengths of the couplings of real 
manganites materials. Obviously, the absence of 3d fer- 
romagnetic phases at half-doping suggests that Jaf is 
substantial in these materials {JafS 2 /t > 0.05). The 
observation of the A-type phase in Pr^S^ / 2 Mn03 with 
distortions that favour the (x 2 — y 2 ) orbitalf£*£ is only 
compatible with g/t > 5.0 (Fig. [2J, a large value that 
is generally corroborated by large distortions^ In addi- 
tion, the fact that many half-doped manganites show a 
(CE) insulator to (ferro) metal transition as a function of 
magnetic- field confirms that g/t can not be much larger 
(since for g/t > 6.3, the ferromagnetic phase is also in- 
sulating). Adopting such values g/t ~ 5 (or E JT ~ t), 
the charge disproportionation in the CE phase is 8 ~ 0.2 
(Fig. |SJ, which is much smaller than Goodenough's ionic 
picture value of 0.5. The inclusion of cooperative JT ef- 
fects is likely to reduce the charge contrast. In addition, 
the specific distribution of the charge contrast amongst 
the Mn and O orbitals depends upon band-structure de- 
tails which we have not considered in this paper 

Secondly, and more importantly, we have studied the 
instabilities of these phases with respect to canting of 
spins and single site defects in their JT distortion pat- 
tern, caused by doping away from x = 0.5 or by the ap- 
plication of a magnetic-field. The consideration of canted 
CE phases allowed us to study how the magnetisation 
changes with applied field. We have found that the dis- 
tortions do not change much in the linear regime of the 
magnetisation up to a threshold field at which there is an 
abrupt change. This seems to be consistent with recent 
experiments^^ A more detailed comparison would be 
interesting in order to extract the strength of the JT cou- 
pling. Regarding the effect of the doping, we have found 
that when electrons are added (with respect to half- 
doping, i.e., for x < 0.5) the transition from an insulating 
CE phase with self-trapped carriers to the ferromagnetic 
metallic phase of the colossal magneto-resistance materi- 
als proceeds via a first-order transition to canted states. 
In contrast, added holes (corresponding to x > 0.5), do 
not favour canted states because of a lack of density of 
states near the top of the valence band. The holes prefer 
to be self-trapped, at least above a threshold g/t. This, 
we believe, is the underlying cause of why the most man- 
ganites tend to remain insulating and favour incommen- 
surate charge orde:»22j2i for x > 0.5, but quickly become 
metallic for x < 0.5. This striking particle-hole asymme- 
try has been recently explored within a Ginzburg-Landau 
framework^!; we believe that our work clarifies the mi- 
croscopic basis for this asymmetry»§i 

We emphasise that the actual numbers we have ob- 
tained, such as the threshold g/t and magnetic-field val- 



ues for the various transitions may not correspond to 
experimental data because the model we have studied 
neglects several important effects such as the coopera- 
tive nature of the JT distortions, the finiteness of the 
Hund's coupling, presence of second neighbour hopping, 
etc. As pointed out in the previous section, for exam- 
ple, the size of the charge gap in the CE phase, which in 
turn determines the threshold field for its transition to 
the ferromagnetic metallic phase, will be reduced when 
these effects are taken into account. 

Finally, and most importantly, in the intermediate JT 
coupling regime which we have argued is the most ap- 
propriate for manganites, we have shown that there is 
an instability of the ferromagnetic phase to formation of 
defects in the JT distortion pattern. We did this by cal- 
culating the energy cost for creating a single site defect 
in the lattice distortion; i.e. reducing (and eventually 
removing) the distortion at that site and promoting a 
quasi-localised electron onto a mobile band. We have 
shown that there are parameter regimes where this ap- 
pears spontaneously. A proliferation of such defects leads 
to a scenario for the ferromagnetic metallic state that is 
completely consistent with, and provides a new justifica- 
tion for, the effective two-fluid (one light and extended, 
the other polaronic and localised) picture proposed re- 
cently to explain the insulator-metal transition in the 
colossal magneto-resistance materials^ 

On the basis of this identification we have suggested 
a new effective Hamiltonian given by eq. which goes 
beyond that of Ref. [33j in that it allows for possible 
orbital and charge ordering effects. We believe that a 
treatment of this Hamiltonian using more sophisticated 
methods, such as the dynamical mean field theory^ will 
eventually lead to a more complete theory of the man- 
ganites, including orbital and charge ordering effects. We 
hope to discuss such work elsewhere. 

Nevertheless, our work suggests that the ferromagnetic 
metallic phase obtained at large magnetic-fields in half- 
doped manganites is similar to the ferro-metallic phases 
found upon hole doping, i.e., for x < 0.5, except, per- 
haps for some remanent orbital and charge order. It is 
an obvious and interesting question as to whether any 
vestige of the orbital and charge order present in the 
CE phase survives metallisation. If it does, the metallic 
phase would also be rather anisotropic, with larger mo- 
bility along the z-axis of the CE phase. But irrespective 
of this, it should have a large fraction of sites which con- 
tinue to be JT distorted accounting for the majority of 
the e g electrons, which remain localised, and only a small 
number of mobile carriers. Interestingly, this resembles 
the phenomenological picture of ref. l36l Our results pro- 
vide a microscopic justification for this picture and can be 
further tested experimentally in a variety of other ways, 
such as measurement of Drude weights in optical conduc- 
tivity, EXAFS and neutron diffraction experiments, for 
instance. 
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APPENDIX A: WAVE-FUNCTIONS OF THE CE 
STATE WITH LATTICE DISTORTIONS 

The unit-cell of the Id zig-zag chains has four inequiv- 
alent sites and there are two orbitals per site (see Fig. 
IB. We l ist below the eight states, which extends Refs. 
|28l29l52fl to the presence of lattice distortions on the 
bridge sites. The corresponding energies are given in Fig. 
ITTland in equations JSJ - (|TU|) . 
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where the notations are t — 4i/3, Ejt — gQ/2 and the 
states are defined by: 



\<i.<\}„. = -^^e^-^l^^a) 



with a = 1, ..,4 the four inequivalent sites of the zig-zag 
chain and a can be any of the orbital states defined in 
eq. flfl. g Q is the component of the wave- vector along the 
chain direction. It takes values in the first Brillouin zone 
[— 7r/2,7r/2]. As for g = 0, the degeneracies at ±7r/2 
come from the translation symmetry combined with a 
mirror plane symmetry. Similarly, the properties \& ga x = 

^aj+7r2 and ^g,,^ = ^q a +TT,4 are consequences of the 
same symmetry. 
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